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ABSTRACT 

Context. Obscured AGN are fundamental to understand the history of Super Massive Black Hole growth and their influence on galaxy 
formation. However, the Compton-thick AGN (Af H >10 24 cm -2 ) population is basically unconstrained, with less than few dozen con- 
firmed Compton-thick AGN found and studied so far. A way to select heavily obscured AGN is to compare the X-ray emission below 
10 keV (which is strongly depressed) with the emission from other bands less affected by the absorption, i.e. the IR band. To this end, 
we have cross-correlated the 2XMM catalogue with the IRAS Point Source catalogue and, by using the X-ray to infrared flux ratio 
and X-ray colors, we selected a well defined sample of Compton-thick AGN candidates at z<0.1. 

Aims. The aim of this work is to confirm the nature and to study one of these local Compton-thick AGN candidates, the nearby 
(z=0.029) Seyfert 2 galaxy IRAS 04507+0358, by constraining the amount of intrinsic absorption (Af H ) and thus the intrinsic lumi- 
nosity. 

Methods. To this end we obtained deep (100 ks) Suzaku observations (A04 call) and performed a joint fit with SWIFT-BAT data. 
We analyzed XMM-Newton, Suzaku and SWIFT-BAT data and we present here the results of this broad-band (0.4-100 keV) spectral 
analysis. 

Results. We found that the broad-band X-ray emission of IRAS 04507+0358 requires a large amount of absorption (larger than 10 24 
cm -2 ) to be well reproduced, thus confirming the Compton-thick nature of this source. In particular, the most probable scenario is 
that of a mildly (Af H ~1.3-1.5xl0 24 cm" 2 , L(2-10 keV)~5-7xl0 43 erg s" 1 ) Compton-thick AGN. 

Key words, galaxies: active - galaxies: individual: IRAS 04507+0358 galaxies: Seyfert - X-rays: galaxies 



1. Introduction 

On the basis of the most accredited X-ray background synthe- 
sis models (Gilli et al. [T997l Treister et al. [20091 Ballantyne 
et al I2006K obscured AGN (Af H >10 22 cm -2 ) dominate the en- 
tire AGN population. For this reason, their space density at dif- 
ferent redshifts is one of the main ingredients in setting the 
evolutionary properties of Super Massive Black Hole (SMBH). 
Unfortunately, the absorption of the obscuring medium (com- 
posed of gas and dust) along the line of sight does not allow 
us to easily detect them and study their nuclear properties in 
the UV/optical bands and in the soft X-ray bands (energy be- 
low few keV). This is particularly true for the most obscured 
sources, the so called Compton-thick AGN, that are hidden by a 
large amount of circum-nuclear obscuring matter along the line 
of sight (column density of A^h>10 24 cm -2 ). 

Compton-thick AGN are generally divided in two main 
classes of sources: mildly Compton-thick with Nu of the order 
of a few times 10 24 cm -2 and heavily Compton-thick with Nu 
above ~10 25 cm -2 . The primary radiation is strongly suppressed 
at low energies in case of mildly Compton-thick AGN, emerging 
only above 10 keV, while in the case of heavily Compton-thick 
AGN the primary radiation is strongly depressed also above 10 
keV, due to the compton down- scattering effect (Matt 1 19961) . In 
both classes of sources, the spectrum below 10 keV is dominated 



by a pure reflection component (i.e. the continuum emission re- 
flected by the putative torus), that is much fainter than the di- 
rect one. For this reason, even for intrinsically luminous objects, 
it is generally difficult to accumulate enough counts below 10 
keV to allow us a reliable spectral analysis and thus to assess 
the nature of these sources. Furthermore, even when the spectral 
analysis is possible, the absorption cut-off is only marginally de- 
tectable below 10 keV (or completely outside the observed en- 
ergy window). This fact does not allow to correctly measure the 
intrinsic absorption and thus to estimate the nuclear luminosity. 
Moreover, as we will show in this paper, in spite of the differ- 
ent values of intrinsic Nu, the shape of Compton-thin AGN with 
5xl0 23 cm" 2 <Af H <10 24 cm" 2 and Compton-thick AGN spectra 
below 10 keV may appear very similar and indistinguishable, 
specially in cases of low counting statistics. This makes even 
more difficult to adequately probe the heavily obscured AGN 
population. Often the presence of Compton thick matter is in- 
ferred through indirect arguments, such as the presence of a 
strong iron emission line at 6.4 keV; for high values of Nu, the 
equivalent width of this line is expected to be high, reaching val- 
ues as a few keV (Matt et al. 1 19961 Murphy & Yaqoob 120091 

To properly set the presence of a Compton-thick source and 
measure the amount of absorption, hard X-ray data above 10 
keV are needed (Comastri et al. 2010 ). These data allow to detect 
the absorption cut-off in the very hard X-ray domain typical of 
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large amount of intrinsic Nu (>10 24 cm" 2 ). For the vast majority 
of Compton-thick AGN without high energy data, only a con- 
servative lower limit on the intrinsic column density and nuclear 
luminosity can be placed. Extreme examples of such limitations 
are NGC 6240 (Vignati et al. [T999l), M rk 231 (Braito et al. l2Q04t 
and Arp 299 (Delia Ceca et al. 120021 where only Beppo-SAX 
PDS (energy range 15-300 keV) observations were able to reveal 
the intrinsic AGN emission, allowing the direct measure of the 
Nu and the intrinsic AGN luminosity. 

A reliable estimate of the number of Compton-thick AGN 
requires first of all the use of efficient methods able to se- 
lect large and possibly complete samples of Compton-thick 
AGN candidates. As a second step, broad-band X-ray spec- 
tra (from a few to hundreds keV) are needed to confirm their 
Compton-thick nature. However, we lack such a statistical sam- 
ple and, more importantly, only a limited fraction of the puta- 
tive Compton thick AGN have been confirmed through broad- 
band X-ray spectroscopy (see Comastri et al. 120041 Delia Ceca 
et all2008al Awaki et al. l200a Teng et al. 2009, Braito et al. 
2009 ). For this reason, indirect arguments are generally used to 
estimate their density in the local Universe (Delia Ceca et al. 
I2008bi) . 

In Severgnini et al. (120091 we found that the (2-10 keV) 
to 24/mi flux ratio vs. X-ray colors can be used as an effi- 
cient technique to select Compton-thick candidates, at least in 
the local Universe. In particular, after having cross-correlated 
the IRAS Point Source Catalog v2.1 (PSC) with the bright end 
(F 2 -io^v>10" 13 erg cm -2 s" 1 ) of the incremental version of the 
2XMM catalogue (Watson et al. 2009), we found that 85% of the 
46 extragalactic sources with F(2-10 keV)/(v24 /im F24 y um)<0.03 
and Hr4^> -0.1 have X-ray properties resembling those of 
a typical Compton-thick AGN. Ten of these are newly discov- 
ered Compton-thick candidates. In order to confirm the nature 
of these objects, we have both cross-correlated our Compton- 
thick AGN candidates with the SWIFT -BAT catalogue (54 
months of observations, Cusumano et al. 120101) and obtained 
Suzaku observations for two of them (IRAS 04507+0358 and 
MCG-03-5 8-007). Here we present the results obtained with the 
broad-band X-ray analysis (Suzaku plus SWIFT-BAT) of IRAS 
04507+0358, the first target observed by Suzaku, while the re- 
sults obtained with the cross-correlation with the SWIFT-BAT 
catalogue will be presented in a forthcoming paper (Severgnini 
et al. in prep.). Throughout this paper we use the cosmological 
parameters Ho =7 1 km s" 1 Mpc -1 , ^=0.7 and Qm=0.3. 



2. IRAS 04507+0358 

This source is optically classified as Seyfert 2 galaxy at z=0.029 
(Strauss et al. 1992). No clear signs of starburst activity have 
been detected (Cid Fernandes et al.l200T). From the IRAS-PSC 
catalogue fluxes, we calculated the infrared (8-1000 /mi) lumi- 
nosity following the prescription reported in Sanders & Mirabel 
(TT9961) and Risaliti et al. (120001) . We found L 7/ ^2.7xl0 44 erg 
s _1 ~7xl0 10 L Q , in agreement with the estimate reported in 
Fraquelli et al. (120031) . 



2.1. XMM-Newton archival data 

The 2-10 keV spectrum used in the analysis described here 
(Observation ID = 0307000401, -11 ksec of good time expo- 
sure in the MOS camera only) was taken from the XMM-Newton 
archive as one of the products of the 2XMM catalogue (Watson 
et al.[2p09 ). The data are grouped with a minimum of 20 counts 
per channel. The X-ray emission of IRAS 04507+0358 (F(2- 
10 keV)~10~ 12 erg cm -2 s" 1 ) could be explained with two for- 
mally acceptable models: a transmission dominated (A^h<10 24 
cm -2 ) and a reflection (A^ H »10 24 cm -2 ) dominated scenario 
(see Fig. 1). In both cases, two thermal components are required 
to fit the softer part of the sp ectrum (belo w 2 keV). W e used 
two mekal models (Mewe et al. [T985l[T986l Liedahl et al. [T995l) 
with the abundances of Wilms et al. ([2000 ), typical of the inter- 
stellar medium. We found CT^O.15^ 4 and £T 2 =0.72+° ^ keV 
(if we use the Compton-thin model) and ^2=0.66^ ^ keV (if 
we use the Compton-thick model). These thermal components 
are most probably associated to star-forming activity. As dis- 
cussed in Section 3, we have estimated a star formation rate 
(SFR) of about 10 M /yr in our source on the basis of the in- 
frared luminosity. For these values of SFR, the presence of a 
multi-temperature components in the X-ray spectra is generally 
observed (e.g. M 82, Persic et al. [20041) . 

As for the data above 2 keV, a clear roll-over is present in 
the spectrum, which, in the Compton-thin hypothesis, could be 
interpreted as a signature of the photoelectric cutoff. We used 
two power-laws with the same photon index fixed to T=1.9 
(Ree ves & Turner [20001 Caccianiga et al. [20041 Page et al. 
2004) to model the scattered and the primary X-ray emission. 
While the scattering component was absorbed only through the 
Galactic column density, we estimated an intrinsic column den- 
sity that absorbs the primary X-ray emission of Nh=(4+2)x10 23 
cm -2 (^ 2 /d.o.f.=41.8/52, left panel of Fig. 1). A similarly good 
fit could be obtained also by replacing the absorbed power- 
law with a reflection component (pexrav model, Magdziarz & 
Zdziarski 1995 ). We call this model the "heavily Compton-thick 
hypothesis" (^ 2 /d.o.f.=53.2/53, right panel of Fig. 1). Even in 
this case, we used the same slope (T=1.9) for both the scattered 
and the reflected component. We fixed the reflection fraction (de- 
fined by the subtending solid angle of the reflector R=Q/27r) 
equal to 1 and the inclination angle to the mean value of 60°. 
We kept the abundances of Wilms et al. ( 2000 ). Independently 
from the assumed model for the underlying continuum, a strong 
narrow Fe line is detected. We fixed the energy line at 6.4 
keV and we found an EW=1.1^ 4 keV (cr=0.2+0.1 keV) in the 
Compton-thin model and EW=2.2+°-jj keV (cr=0.3+°- 4 keV) in 
the heavily Compton-thick hypothesis. 

Apart from the strong iron line that could be considered as 
an indirect hint for the presence of a Compton-thick source, 
this analysis shows that the XMM-Newton data alone can not 
distinguish between a Compton-thin or Compton-thick scenar- 
ios. With the aim of assessing the nature of this object and 
estimating its intrinsic properties (e.g. Nu and luminosity), we 
have obtained 100 ks of Suzaku observation in the A04 call. We 
have combined here this Suzaku data with SWIFT-BAT spec- 
trum accumulated during the first 54-months of observations 
(Cusumano et. 12010b . 



1 HR4 is defined using the two following bands: 2-4.5 keV and 4.5- 
12 keV: HR4= g|g|^g|f§E|i^ , where CTS are the vignetting 
corrected count rates in the energy ranges reported in bracket. See 
Watson et al. l2009l for details. 



2.2. Suzaku data 

IRAS 04507+0358 was observed by the Japanese X-ray satel- 
lite Suzaku (Mitsuda et al. 2007 ) for 100 ks at the beginning of 
September 2009. At the time of these observations only three 
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Fig. 1. XMM-Newton (MOS1 - black data points in the electronic version and MOS2 - red data points in the electronic version) 
unfolded spectrum of IRAS 04507+0358 fitted with a Compton-thin (left panel) and a Compton-thick model (right panel). In the 
lower panels the ratio between data and best-fit models are shown. In the Compton-thin hypothesis, the data are best reproduced 
by two thermal emission components (a, b) plus a power-law absorbed only by the Galactic column density (scattering component, 
c) and a power-law absorbed by a high column density (Nh=(4±2)x10 23 cm -2 ) absorber (e). The two power-laws have the same T 
fixed to 1.9. A prominent iron line is also present (d). A similarly good fit could be obtained by using a cold reflection component 
instead of the absorbed power-law (e component in the left panel). 



of the four X-ray Imaging Spectrometers (XIS, Koyama et al. 
120071) were working: two front-illuminated (XIS0 and XIS3) 
and one back-illuminated (XIS1) CCDs. For our analysis we 
used the 0.4-8 keV and 0.4-10 keV data obtained by XIS1 and 
XIS0+XIS3 respectively, combined with the HXD-PIN data. 
This latter is a non imaging hard X-ray detector (Takahashi et al. 
120071) covering the 12-70 keV energy band. IRAS 04507+0358 
was placed at the HXD nominal pointing. Cleaned event files 
were filtered with the standard screening The effective expo- 
sure time after data cleaning are 83.7 ks for each of the XIS 
source's and 77.6 ks for the HXD-PIN. 

XIS spectra: were extracted from a circular region of 2.3 
arcmin of radius centered on the source. Background spectra 
were extracted from two circular regions with the same radius of 

2 The screening procedure filter all events within the South Atlantic 
Anomaly (SAA) as well as with an Earth elevation angle (ELV) <5° and 
Earth day-time elevation angles (DYE_ELV) less than 20°. Furthermore 
also data within 256s of the SAA were excluded from the XIS and 
within 500s of the SAA for the HXD. Cut-off rigidity (COR) criteria 
of >8GV for the HXD data and >6GV for the XIS were used. 



the source region but offset from the source and the calibration 
sources. The XIS response (rmfs) and ancillary response (arfs) 
files were produced using the latest calibration files available and 
the ftools tasks xisrmfgen and xissimarfgen, respectively. The 
net count rates observed with the three XIS in the 0.4-10 keV 
band are 0.034+0.001 cts/s (XIS0), 0.048+0.013 cts/s (XIS1), 
0.040+0.001 cts/s (XIS3). 

The spectra from the two front-illuminated CCD were then 
combined, while the back-illuminated CCD spectrum was kept 
separate and then fitted simultaneously. The net XIS source spec- 
tra were then binned in order to have a minimum signal-to-noise 
ratio (S/N) ^4 in each energy bin. 

HXD-PIN spectrum: At the time of writing, two instru- 
mental background (called non-X-ray background, NXB) files 
have been released, the quick and the tuned one. We tested both 
of them and we found that the tuned background count rate in 
the 15-70 keV is 5% higher than the count rate of the quick 
one in the same energy range. To estimate which is the event 
file that provides the most reliable estimate of the real NXB, we 
compared them with the data taken during periods of Earth oc- 
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cultation in the same energy range, 15-70 keV. Since the Earth 
is known to be dark in hard X-rays and during the Earth oc- 
cultation we do not measure the emission of the source and of 
the cosmic X-ray background, the occulted data should give a 
good representation of the actual NXB rate. Thus, after having 
corrected the Earth occultation data for the dead-time and taken 
into account only the events with an Earth elevation angle lower 
than -5°, we extracted the Earth occultation spectrum and light- 
curve in the 15-70 keV an 15-50 keV, respectively. 

The comparison between the occulted data and the quick 
(bkgA) and tuned (bkgD) backgrounds is shown in Fig. 2. The 
tuned background level is systematically higher than the quick 
background level and than the Earth occultation level. This dif- 
ference becomes larger at E>35 keV. The light-curve compar- 
ison is shown in Fig. 3. The tuned background spectrum and 
light-curve are nearly 7% above the data taken during the Earth 
occultation. 

As a final check, we compared Suzaku flux with Swift-BAT 
data. We found that with the quick background, the Suzaku flux 
is close to the Swift one, i.e. F( 15-70 k&V)suzaku= 1.9x10" 11 erg 
cm -2 s" 1 andF(15-70keV) fiAr =1.8xlO -11 erg cm -2 s" 1 (assum- 
ing the same model for the two instruments). We thus decided 
to use the quick background file and combined it with the cos- 
mic X-ray background, the latter was parametrized by us using 
the prescription of Boldt STWJl and Gruber et al. ST999i . The 
background-corrected count rate in the 15-70 keV is 0.033 cts/s 
(-8.3% above the background level). For the spectral analysis, 
we rebinned the HXD-PIN spectrum in order to have a signal- 
to-noise ratio of 4 in each energy bin. 

2.3. SWIFT-BAT data 

The SWIFT-BAT data were processed using the BatJmager 
software (Segreto et al. [20T0b . IRAS 04507+0358 was detected 
with a S/N~12 in the 14-150 keV energy range. In order to pro- 
duce the spectrum of IRAS 04507+0358 integrated over the 54 
months of survey data, we extracted the rates and the relevant 
errors from the pixel corresponding to the best BAT position in 
the all-sky maps produced in eight energy bands (14-20 keV, 20- 
24 keV, 24-35 keV, 35-45 keV, 45-60 keV, 60-75 keV, 75-100 
keV, 100-150 keV). The spectrum was analyzed using the BAT 
spectral redistribution matrix distributed with the BAT CALDB. 
Further details can be found in Cusumano et al. ([2009 ). 

2.4. Broad band X-ray spectral fitting 

The XIS, the HXD-PIN and the BAT spectra have been fitted 
simultaneously covering a wide energy range (0.4-100 keV) 
using Xspec version 12.5.0. We tied together XIS, HXD and 
BAT parameters. We left free to vary the BAT/XIS cross- 
normalization, while for the HXD/XIS instruments we assumed 
a cross-calibration of 1.18 (Manabu et al. 2007; Maeda et al. 
200^}. In all the models described here we have used the abun- 
dance of Wilms et al. ( 2000 ) and the Galactic hydrogen column 
density along the line of sight (from Dickey & Lockman ll990K 
Af H (Galactic)=6.7xl0 20 cm" 2 . 

Since the aim of this work is to explore the physical proper- 
ties of the nuclear regions of IRAS 04507+0358 (e.g. the N u 
of the obscuring matter and the de-absorbed X-ray luminos- 
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Fig. 2. Comparison between the quick (bkgA, red triangles) and 
the tuned (bkgD, green crosses) backgrounds with the Earth oc- 
cultation spectra (black squares). The BkgD spectrum is above 
both the bkgA and the Earth data. The tuned background clearly 
overpredicts the real NXB, in particular above 35 keV. 
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Fig. 3. Comparison between light curves extracted in the 15-50 
keV energy range. The Earth light curves are corrected for the 
dead time. In the upper panel we show the difference between 
the Earth and the tuned background (bkgD) light curves, while 
in the lower panel we show the difference between the Earth and 
the quick background (bkgA) light curves. 



ity of the AGN), a detailed analysis of the spectrum below 2 
keV will be not discussed here. Starting from the results ob- 
tained with XMM-Newton data, we have fitted this part using 
two thermal components. We left the temperatures free to vary 
and we found that the two temperatures are in good agreement 
with those found with XMM-Newton data: CTi=0.14^ and 

*T 2 =0.73^;gJ. 



Severgnini et al.: Suzaku and SWIFT-BAT observations of a newly discovered Compton-thick AGN 



Unfolded spectrum and Compton-thin model 
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Fig. 4. Suzaku XIS-FI (black data points in the electronic 
version), XIS-BI (red in the electronic version), HXD-PIN 
(green in the electronic version) and SWIFT-BAT (blue in the 
electronic version) data of IRAS 04507+0358. Upper panel: 
Unfolded Compton-thin model. The model consists of: two ther- 
mal emission components (a and b components), a scattered 
component (c), a Fe emission line at -6.4 keV (d) and an ab- 
sorbed power-law with a Nu fixed to 4xl0 23 cm -2 (e). Lower 
panel: Ratio between data and the Compton-thin best-fit model. 



As for the XMM-Newton spectra, to fit the data above 2 
keV, we used a simple basic model for obscured AGN, i.e. two 
power-laws with the same slope. One of the two power-laws 
is absorbed only by Galactic column density and represents 
the scattered component; the second power-law represents the 
primary X-ray emission seen through the cold absorber (the 
putative torus). To reproduce this absorbed power-law, we used 
the model by Yaqoob 119971 (plcabs in Xspec), which partially 
takes into account for the Compton down- scattering. The Fe 
emission line was modeled with a Gaussian profile. 

Compton-thin hypothesis: As a first step, we tested the 
Compton-thin hypothesis by fixing the slope of the power-laws 
to r=1.9 and the column density to A^h=4x10 23 cm -2 , i.e. 
the values found with the XMM-Newton data (see Fig. 4, 
;^ 2 /d.o.f.=527.5/194). This model provides an overall good rep- 
resentation of the XIS data, but it under-predicts the HXD-PIN 




Energy (keV) 

Fig. 5. Residuals, with respect to a single power-law, of the 
XIS data/model (XIS-FI - black data points in the electronic 
version and XIS-BI-red points in the electronic version) of 
IRAS 04507+0358 at energies close to the Fe band (observed 
frame). No iron line is included in the model. 



and the BAT data. We left the photon index and the column den- 



— sity free to vary and we found r= 1.8^9 an d Nu 



=5.5^; 2 xl0 23 

cm -2 (x 2 /d.o.f. =334.95/1 93); however also this model does not 
provide a good representation of the emission above 10 keV. 
Our first conclusion is that the Compton-thin hypothesis is 
clearly ruled out by the HXD-PIN and BAT data. 

Mildly Compton-thick hypothesis: To account for the 
prominent Fe emission line in the spectrum and for the excess 
detected above 10 keV, we added to the model a component, 
which represents the emission reflected from neutral material 
into our line of sight. In particular, we added the pexrav model 
leaving open the possibility that the reflected emission could be 
partially obscured by the torus itself. We fixed the inclination 
angle to its default value (i = 60°) and the normalization equal 
to that of the intrinsic powerlaw. With the addition of this new 
component, that we found to be absorbed (Nh=(3.3±0.2)x10 23 
cm -2 ), the model provides a relatively good fit to the 0.4-100 
keV spectrum 0f 2 /d.o.f. = 265.5/190). We found a BAT/XIS 
cross-normalization of l.l^or The primary X-ray continuum 
has an intrinsic slope of T=2.4±0.5 and is absorbed by an 
Af H =(1.45±0.10)xl0 24 cm -2 . We found that the scattering frac- 
tion is less than 1%, and the reflection fraction is R=0.7^. The 
slightly steep photon index, which is also not well constrained 
(ranging from 1.9 to 2.9), is not so striking if we take into ac- 
count the complexity of the model. 

The Suzaku data confirm the presence of a Fe Ka line at 
6.37+0.01 keV (<x<50 eV), but with a lower equivalent width 
(EW=450^q° eV) with respect to the value measured with the 
XMM-Newton data, although consistent considering the errors. 
In order to investigate the presence of other possible features, we 
inspected the residuals around 6.4 keV with respect to a single 
power-law model (see Fig. 5); there are clear residuals, both red 
and blue-wards the Ka line. Some excesses in the data/model 
ratio are present around 5.8 keV (observer's frame) and between 
6.5 and 7 keV, while at E> 7 keV there is a clear drop. While 
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the first excess is most probably due to a bad subtraction of the 
emission lines present in the calibration source (i.e. Mn Kal 
at 5.899 keV and Mn Ka2 at 5.888 keV), the drop at energies 
larger than 7 keV is most likely due to a high level of the XIS 
background. Finally, the excess between 6.5 and 7 keV could 
suggest the presence of a second He- or H-like Fe emission 
line. In particular, we found that a second narrow (<x<50 eV) 
line centered at -6.92 keV rest-frame (EW-100 eV) could be 
actually present (see Fig. 6) even if it is not statistically required 
(^ 2 /d.o.f.=255.9/188). This second line could be associated both 
to the Fe XXVI Lya emission line (rest-frame energy E=6.96 
keV) or to the Fe KJ3 emission line (rest- frame energy E=7.058 
keV) or to a combination of these two lines. However, we es- 
timated the flux of the 6.92 keV line and found that it is about 
15% of the flux of the Ka line, as it is expected for the Fe KJ3 
emission line (see e.g. Leahy & Creighton ll993b . 

As reported in Table 1, the best-fit values obtained with 
both models (one and two emission lines) are consistent within 
the errors. We refer to these models as mildly Compton-thick 
AGN models. In Table 1, the lines at ~6.4 and ~6.9 keV are 
marked as Ei and E2 respectively. The 2-10 keV flux, corrected 
for the Galactic absorption along the line of sight, is F(2-10 
keV)=2xlO" 12 erg cm" 2 s" 1 and, once corrected also for the 
amount of intrinsic absorption, the deabsorbed 2-10 keV lumi- 
nosity of the AGN is 7xl0 43 erg s" 1 . 

To further test the above models, we modeled the two com- 
ponents associated with cold reflection (the 6.4 keV Gaussian 
emission and the Compton- scattered continuum from neutral 
material) using a single and self-consistent model: reflionx (Ross 
& Fabian 2005 ). In this model, an optically thick disc is illumi- 
nated by a power-law, producing florescence lines and contin- 
uum emission. The parameters include the Fe abundance, ion- 
ization parameter £ (defined as £ = 47rFtot/n#, where Ftot is 
the total illuminating flux and n H is the density of the reflec- 
tor) and the incident power-law photon index T. Being the ob- 
served Fe iron Ka emission line unresolved in the Suzaku spec- 
trum, no additional velocity broadening was applied to the re- 
flected spectrum. We found that the data are well reproduced by 
this model (^ 2 /d.o.f.=226. 3/190), confirming the physical con- 
sistence between the measured reflected continuum and the line 
intensity. In this case the reflected component is absorbed by 
A^h=(2.7±0.8)x10 23 cm -2 , while the transmitted component is 
absorbed by a larger column density A^h=(1.3±0.3)x10 24 cm -2 . 
We found a value for the photon index that is consistent with the 
previous ones (r=2.5±0.5) and a ionization parameter £=29^ 3 
erg cm s" 1 (the lowest value allowed by the model is £=10 
erg cm s" 1 ), in agreement with iron atoms typically in a low- 
ionization state corresponding to Fe I-XVII. Best fit parameters 
are reported in Table 1 (mildly Compton-thick AGN models - 
reflionx). 

Finally, we checked the self-consistence of the mildly 
Compton-thick hypothesis by using also the recent paper pub- 
lished by Murphy & Yaqoob (12009b . They calculated Green's 
functions that may be used to produce spectral- fitting routines to 
model the putative neutral toroidal X-ray reprocessor in AGNs 
for an arbitrary input spectrum. In their calculation the repro- 
cessed continuum and fluorescent line emission due to Fe Ka, Fe 
Kfi and Ni Ka are treated self-consistently. On the basis of the 
EW obtained for the Fe Ka emission line (~450 eV) by our anal- 
ysis and by assuming an inclination angle between the observer's 
line of sight and the symmetry axis of the torus larger than 
60°, the model predicts an Af H -1.5-3xl0 24 cm -2 . The mildly 
Compton-thick hypothesis is thus fully supported also by this 



recent model proposed in the case of neutral toroidal X-ray re- 
processor in AGNs. 

Heavily Compton-thick hypothesis: For completeness, we 
tested also the heavily (A^ H >10 25 cm" 2 ) Compton-thick AGN 
hypothesis. In this case, the emission is completely dominated 
by the scattered component at low energies and only by the re- 
flected component (reflionx or pexrav model) at high energies. 
By using disc-reflection models (reflionx or pexrav models), we 
found that this scenario is statistically acceptable. As an exam- 
ple, we report in Fig. 7 and Table 1 the results found with the 
reflionx model Or 2 /d.o.f. = 218.9/189, heavily Compton-thick 
AGN models - reflionx in Table 1). However, an intrinsic column 
density higher than 10 25 cm" 2 is not supported by the model of 
Murphy & Yaqoob ([20091) for Fe Ka EW of order of 450 eV. 
Moreover, on the basis of this model, the observed luminosity 
should be just less than few percents of the intrinsic one, by im- 
plying an intrinsic luminosity of L(2-10 keV)>10 45 erg s" 1 . This 
value of X-ray luminosity exceeds the infrared luminosity mea- 
sured for this source (L/#=2.7xl0 44 erg s" 1 ), which, under the 
assumption that most of the optical and ultraviolet radiation is 
absorbed by a dusty torus surrounding the nuclear source, is a 
good proxy of the bolometric luminosity. Therefore, on the ba- 
sis of the Murphy & Yaqoob (2009 ) model, this scenario is not 
acceptable from the physical point of view. 



3. Discussion and Conclusion 

IRAS 04507+0358 was proposed as a new Compton-thick 
AGN candidate through the diagnostic diagram described in 
Severgnini et al. (2009 ). Here, by studying its broad band X- 
ray spectral properties (0.4-100 keV), we confirm the Compton- 
thick nature of this AGN. 

We showed in this paper that with the XMM-Newton and 
XIS data below 10 keV alone, we could not distinguish be- 
tween Compton-thin or Compton-thick scenarios. However, the 
Compton-thin model clearly under-predicts higher energy data. 
In particular, we found that the most favorite scenario is that 
of a mildly Compton-thick AGN with N H =1.3-1.5xl0 24 cm" 2 
and L(2-10 keV)=5-7xl0 43 erg s" 1 . This luminosity is in full 
agreement with the observed infrared luminosity and with that 
predicted for the torus by Fraquelli et al. (2003 ). These authors, 
after having estimated the rate of the ionizing photons emitted 
by the AGN and assuming an opening angle of 30° for the ion- 
ization cones, predicted an infrared luminosity for the torus of 
L/#=5.3xl0 43 erg s" 1 . This value is about 20% of the infrared 
luminosity (L/#^2.7xl0 44 erg s" 1 , see also Sect. 2). The remain- 
ing fraction of the infrared luminosity (2.2x1 44 erg s" 1 ) is most 
probably due to star- formation activity and it is roughly in agree- 
ment with the soft X-ray luminosity of the thermal components 
of our source (L(0.5-2 keV)=2.5xl0 41 erg s" 1 ) considering the 
L(soft-X)-L(FIR) relation presented in Persic et al. (120041) . To 
estimate the relevant SFR of our source, we adopted the follow- 
ing Kennicutt et al. (fT998t relation: SFR=(IW5.8xl0 9 L Q ) M Q 
yr" 1 and we considered only the infrared luminosity produced by 
star-formation activity (~2.2xl0 44 erg s" 1 ). We found a SFR of 
about 10 Mo/yr. 

As for the heavily Compton-thick hypothesis, by using disc- 
reflection models it is possible to well reproduce the broad-band 
spectrum of IRAS 04507+0358 and to obtain a statistically ac- 
ceptable fit of the data. In this case, in order not to exceed the 
bolometric luminosity of the source, the reflection component 
should be less than 10% of the intrinsic one, which is not sup- 
ported by the Murphy & Yaqoob (2009 ) model. 
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Table 1. Best-fit values of the mildly and heavily Compton-thick AGN models 



Mildly Compton-thick AGN models 


Model 


r 


N H 
[10 22 cm" 2 ] 


Refl. 
frac. 


Scatt. 
frac. 


Ei 
keV 


cr 
eV 


E 2 
keV 


eV 


*7d.o.f. 


L(2-10keV) 
10 42 erg s" 1 


One line 
Two lines 


2.4+0.5 
2.4+0.5 


145+10 
1463i 


7 +04 
-Q-5 
7 +0 - 4 


<1% 
<1% 


6.37+0.01 
6.37+0.01 


<50 
<50 


6 92 +0 - 50 


<50 


265.5/190 
255.9/188 


70 
70 


Model 


r 


Nu 
[10 22 cm" 2 ] 


erg cm s 1 


E 2 
keV 


o~ 

eV 


^7d.o.f. 


L(2-10keV) 
10 42 erg s" 1 








Reflionx 


2.5+0.5 


130+30 


29!f 


6 91 +012 


<50 


226.3/190 


50 








Heavily Compton-thick AGN model 


Model 


r 


[10 22 cm" 2 ] 


? 

erg cm s 1 


E 2 
keV 


cr 
eV 


X z /d.o.f. 


L(2-10keV) 
10 42 erg s" 1 








Reflionx 


2.2+0.5 


>100(Q 


66 + _)° 


6 93+ - 07 


<50 


218.9/189 


>10 3 









a This value of Nu is not an output of the fit, but it is an assumption that we made in this model. 



We note that this source could not be found as a possi- 
ble Compton-thick AGN by using only the data below 10 keV, 
nor by using the two-dimensional diagnostic tool for reflection- 
dominated Compton-thick object proposed by Bassani et al. 
(1999). This latter diagram shows that Compton-thick AGN 
should be characterized by high Ka iron line equivalent width 
(EW>300 eV) and by a 2-10 keV flux normalized to the [OIII] 
optical-line flux (T parameter) typically lower than 0.5. In this 
estimate, the [OIII] optical-line flux should be corrected for 
Galactic and intrinsic extinction. From the X-ray analysis per- 
formed on IRAS 04507+0358, we derived for the Fe Ka line an 
EW of -450 eV. In order to estimate the T parameter, we con- 
sidered the [OIII] line flux (already corrected for the Galactic 
extinction) reported by Cid Fernandes et al. (2001). After the 
correction for the intrinsic extinction^] we found T=2, larger 
than the typical value adopted to select Compton-thick candi- 
dates («1). This implies that by using 2-10 keV data alone 
and/or diagnostic diagrams based on the Fe Ka line EW vs. 2- 
10 keV to [OIII] flux ratio, a significant fraction of mildly or 
heavily Compton-thick AGN could be missed. Whereas some 
examples of Compton-thick AGN with low value of Fe Ka 
line equivalent width are indeed already present in the literature 
(e.g. Awaki et al. [20001 Ueda et al. 120071 Braito et al. 1 20091). 
IRAS 04507+0358 escapes the diagnostic diagram due to its 
high value of T. In the case of mildly Compton-thick AGN, 
this is most probably due to the large amount of light reflected 
from the torus (i.e. about 35% of the light emitted by the central 
source) that increases the continuum in the 2-10 keV range. The 
Compton-thick AGN missed using the Bassani et al. (11999b cri- 
teria could be selected by combining infrared information with 
2-10 keV data and with very hard X-ray follow-ups, as shown 
by the study presented here. Moreover, our analysis, besides 
providing a further evidence of the potentiality of the F(2-10 
keV)/F24yum vs. HR diagnostic diagram, also shows the impor- 



4 For the intrinsic extinction, we applied the following formula: 
Fcor[OIII] = Fobs[OIII][(HQf/H^)obs/(HQf/H^) ] 2 - 94 (Bassani et al. 
[1999}, where Fcor[OIII] is the extinction-corrected flux of [OIII]5007, 
and Fobs [OIII] the observed flux of [OIII] 5007. We used and observed 
line ratio of Hof/H/?=5.5 (Cid Fernandes et al. 2001) and an intrinsic 
B aimer decrement of (Ha/H/?) =3.0. 



tance of using a wide X-ray spectral coverage in order to con- 
strain the intrinsic column density in this type of sources even in 
the case of mildly Compton-thick AGN. 

With the aim of establishing the nature of other Compton- 
tick AGN candidates found through our diagram, we have re- 
cently obtained in the Suzaku A05 call 100 ksec of observation 
for MCG-03-58-007, while the broad-band X-ray analysis for 
other twelve sources with a SWIFT-BAT counterpart is already 
ongoing and will be the subject of a forthcoming paper. The fi- 
nal aim is to better constrain and define the Compton-thick AGN 
population and hence to better estimate their space density. 
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Fig. 6. Suzaku XIS-FI (black data points in the electronic 
version), XIS-BI (red in the electronic version), HXD-PIN 
(green in the electronic version) and SWIFT-BAT (blue in the 
electronic version) data of IRAS 04507+0358. Upper panel: 
Unfolded mildly Compton-thick model. The spectral compo- 
nents are: a) and b) two thermal emission components; c) scat- 
tered AGN component; d) narrow Gaussian line at 6.37 keV; 
e) absorbed power-law AGN component; f) pure reflection 
AGN component; g) narrow Gaussian line at -6.9 keV. Lower 
panel: Ratio between data and the mildly Compton-thick best- 
fit model. 
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Fig. 7. Suzaku XIS-FI (black data points in the electronic ver- 
sion), XIS-BI (red in the electronic version), HXD-PIN (green 
in the electronic version) and SWIFT-BAT (blue in the elec- 
tronic version) data of IRAS 04507+0358. Upper panel: The 
unfolded heavily Compton-thick model (reflionx, see Table 1) 
is overplotted on the data. The spectral components are: a) and 
b) two thermal emission components; c) scattered AGN com- 
ponent; d) narrow Gaussian line at ~6.9 keV; e) pure reflection 
AGN component (reflionx). Lower panel: Ratio between data 
and the heavily Compton-thick best-fit model. 
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